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Abstract 

The double differential multiplicities and rapidity distributions for A hyperon pro- 
duction in central Au + Au interactions at AGS in the range of rapidities from 
1.7 to 3.2 and the range of transverse kinetic energies from 0.0 to 0.7 GeV are 
parametrized in terms of the the Blast Wave approximation. The longitudinal and 
transverse radial expansion parameters and the mean temperature of A hyperons 
after the freeze-out of the nuclear matter are presented. The predictions of the 
RQMD model with and without mean field potentials are compared to our data. 
Both variants of RQMD are parameterized in terms of the Blast Wave model and 
the results of such parameterizations are compared to the experimental ones. It is 
found that inclusion of the mean field potentials in RQMD is essential to account 
for the strong expansion observed in the data. 



1 Introduction 

The collective nuclear flow has been a subject of extensive theoretical and 
experimental investigations for the last few years. It was pointed out that, 
the collective flow directly probes the equation of state of nuclear matter and 
thus can provide one with information about the possible phase transition 
[1,2]. Before a conclusion about the phase transition can be made, the nuclear 
equation of state should be known for the baryon densities far above normal, 
but where the matter is still in the hadronic phase. However the theoretical 
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statements for the compressibility of hadronic phase or for the stiffness of 
the equation of state contradict each other [1,3]. Measurements of collective 
nuclear flow may specify the limits of possible compressibilities of hadronic 
matter. 

One usually distinguishes between longitudinal flow, transverse directed flow, 
elliptic flow and transverse radial expansion, which are interrelated forms of 
a global picture [3,4]. The directed flow comprises a correlation between the 
impact parameter vector and the directions of produced particles. For quasi- 
central heavy-ion interactions this correlation vanishes and one encounters 
the azimuthal symmetry of flow with respect to the impact parameter: the 
transverse radial expansion superposed on the longitudinal expansion. The 
spectra of particles produced in the heavy-ion interactions represent the con- 
volution of stochastic thermal and collective macroscopic motion. The shapes 
of one-particle multiplicities carry part of the information on the collective 
velocity profile together with information on the stochastic fluctuations of the 
particle velocities. Even though it is impossible to deconvolve the spectra pre- 
cisely, one can obtain a semi-quantitative picture of nuclear matter expansion 
with certain assumptions about the flow velocity profile. Studying the spec- 
tra of A hyperons has several advantages over studying the spectra of pions 
and protons. First of all, in central heavy-ion collisions the majority of these 
particles are produced via multi-step interactions, due to the suppression of 
the strangeness production in nucleon-nucleon collisions. Thus, the A hyperon 
spectra essentially miss a non-hydro dynamical component, which is present 
in the spectra of protons. Further, unhke the tt mesons the influence of reso- 
nances on the shapes of the A spectrum is negligible. And finally, being heavy, 
the A hyperons acquire significant momentum due to the collective component 
of their velocities. 

In this paper we will show that the double differential multiplicities for the 
A production in central Au + Au collisions at AGS, experiment E891, can 
be described in terms of an aspherical blast wave with a single temperature. 
We will extract the parameters of the blast wave from our data and from 
the predictions of two versions of the RQMD model (Relativistic Quantum 
Molecular Dynamics [3]): the cascade version and the version with the mean 
fields turned on. 



2 Experimental method 

We analyze the experimental data of the E891 collaboration on the A hyperon 
production in central Au + Au collisions at the BNL AGS accelerator [10]. The 
plan view of the E891 apparatus is shown in Figure 1. The 11.6 x A GeV/c 
Au beam was incident on a 0.025 cm Au target (7.5% radiation lengths). The 
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charged particles were tracked in three Time Projection Chambers (TPC) in 
the 10 kG field of the Multiple Particle Spectrometer (MPS) magnet. The 
field volume of the MPS magnet is 4.60 meters long xl.80 meters wide xl.20 
meters high. Downstream of the TPC's were four Multiwire Drift Chambers 
used for track calibration. The Au target was 80 cm upstream of the active 
volume of the first TPC. This distance was chosen to maintain adequate two 
track separation to reduce track merging in the front TPC module. The beam 
was incident on the center of the first TPC. This arrangement covered the 
forward hemisphere in the NN center of mass system for rapidities > 1.7. 
The beam intensity varied from 500 to 2000 Au ions per second. 

The drift volume of each TPC module is 60 cm high x65 cm wide x47 cm 
long. The drift volume was filled with a low-diffusion gas mixture (79% Argon, 
16% Isobutane and 5% Dimethoxymethane) known to be stable at high gain. 
The TPC's were operated at 20 kV. The high voltage was apphed parallel to 
the magnetic field. The three dimensional coordinates of the hits in the drift 
volume of each TPC were measured by 12 rows (spaced 3.8 cm apart in the 
^-direction) of 1 cm long anode wires with 256 wires in each row, (2.54 mm 
spacing in the x direction). The drift time provides the vertical coordinate, 
y, measured in 1024 time bins and the wire number gives the two horizontal 
coordinates. The z axis points along the beam, the x axis is horizontal and 
points towards the bend direction of positive particles, and the y axis points 
up with the magnetic field directed along the y axis. Electrons produced by 
the charged particles in the drift region must pass through a gating grid and a 
wire cathode in order to reach the anode wires. A second, solid copper cathode, 
below the anodes helps to achieve the high gain (Figure 2). The use of short 
anode wires precludes making dE/dx measurements. Further details about 
the TPC construction are given in [5]. 

Beam intensity limitations were imposed by positive ion buildup in the drift 
region and in the amplification region of the TPC. ([7,8]). For the ions created 
in the amplification region, the gating grid was used to prevent the transfer 
of drifting electrons from the drift region to the amplification region and that 
of positive ions in the reverse direction except for selected events. Gating 
was achieved by applying different potentials to adjacent wires, as shown in 
Figure 2 collecting the electrons on the more positive wires. When the trigger 
system indicates that an event of interest occurred the gate wires were shorted 
together for 40 ms, a little longer than the maximum drift time of electrons 
(the drift velocity of e~ is 1.6 cm/ms), allowing ionization from the tracks to 
reach the amplification region. Because of their slower drift velocity, positive 
ions from the amplification at the anode do not reach the gating grid before 
the gate is closed, further improving the gating of the TPC. 

The electronics were mounted directly on the chamber. The circuit for each 
channel, containing an amplifier-shaper-comparator and memory, is shown 
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Fig. 1. The E891 plan view. Dimensions are out of scale. 
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Fig. 2. Schematic arrangement of electrodes in the detector end cap. 

in Figure 3. Sixteen channels were assembled in a hybrid package (LeCroy 
HTD 161S,M) with inputs spaced 2.54 mm apart. Shaping time constants 
were selected to match the anode waveform. An integration time constant 
of 50 ns and a differentiation time constant of 200 ns were chosen in order 
to minimize noise and maximize two track separation. Each hybrid contains 
memory capable of recording 1024 time samples from each channel. To reduce 
local noise pickup, all logic levels are ECL. The clock and address lines are 
terminated differential ECL. The on-chamber performance of the hybrid allows 
input sensitivity of less than 1 ^A. One MASTER and seven SLAVE hybrid 
circuits were assembled on one printed circuit board which is sufficient for 
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Fig. 3. Scheme of individual electronics channel. 

128 anode readouts and is externally accessed by one clock line, one data line 
and 14 control lines. The Master hybrid had the control logic necessary to 
read out itself plus the seven slave hybrids. The time samples were obtained 
at 32 MHz and the data is sequentially read from each channel's memory at 
16 MHz. This serial string of data bits is encoded into wire number, drift 
time and cluster size by FASTBUS modules in which all rows are encoded 
simultaneously. Consequently the entire readout is completed in 10 ms (dead 
time). More details on the FASTBUS system are given in [9] 

The position of the incident beam was measured by two proportional wire 
chambers upstream of the MPS magnet. The trigger selected centrally enriched 
events for data recording. A plastic Cherenkov counter CI — C2 was used 
to select Au ions. A halo veto scintillation counter 5*2 was used to reject 
events with interactions upstream of the target. The beam logic was BEAM — 
SI- S2- CI ■ S3. In order to enrich the trigger for central events, a multiplicity 
coimter 5*4 was used to select events in which the beam particle interacted with 
the target and the counter S6 vetoed fragments with Z > 6. Si was a lOx 10 m 
scintillation paddle placed just above the beam. The logic for a trigger was 
TRIGGER = BEAM ■ S4:-SE. The pulse height in S6 was recorded for 
every event so that tighter cuts can be applied in the software. Multiplicites 
of reconstructed tracks often exceed 200, with an estimated reconstruction 
efficiency > 90% for long tracks. 

The drift chambers downstream of the TPC's were used for track calibration 
of the TPC's. A target was placed in the beam line upstream of the apparatus 
to generate a flux of fast parallel tracks fllling the TPC volume. These tracks 
were measured in the TPC's and drift chambers for field on and off conditions. 
Extrapolating the drift chamber tracks to the TPC hits under these conditions 
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were used to determine the alignment of the TPC's (field off) and the E x B 
corrections (field on) due to electric field inhomogeneities at the edges of the 
TPC's. 



2.1 Data analysis 

The data and analysis of this experiment were published previously in [10]. 
Events with multiplicities of charged tracks ranging up to 280 were recon- 
structed, using a three-dimensional tracking program. The total charged track 
multiplicity at the interaction vertex as reconstructed by the tracking program 
was used as a measure of centrality. We found from Monte-Carlo studies that 
this correlates reasonably well with the impact parameter. For the results 
shown in [10] we selected events with total charged track multiphcities in the 
acceptance of >220, corresponding to a measured cross section of 270 mb 
(or ~ 4.4% of the geometric Au + Au cross section). This selection criterion 
resulted in a sample of 14,114 central events. 

Tracks which missed the interaction point by more than 7 mm were taken 
to be candidates for a decay vertex. To reduce combinatorial background, we 
required the decay point to be more than 15 cm downstream of the interaction 
point, the reconstructed momentum vector to point back to the production 
point within 3.5 mm (3 a cut), and the decay plane to be out of the mag- 
netic field bend plane. In addition, in order to select tracks which had good 
momentum resolution, we required the sagitta of the measured tracks to be 
> 0.67 cm. The data was recorded with a beam rate of ~ lO^Au ions/sec 
(~ 6 X 10^ minimum ionizing). This beam rate resulted in noticeable space 
charge distortions in the TPC modules in the beam region [7,8], requiring 
us to remove from the data sample A hyperons decaying with the proton in 
this region. Specifically, the A hyperon candidates with protons having an az- 
imuthal angle within the range, —52° < ip < 149°, were removed, where the 
azimuthal angle, (p, is measured with respect to y axis. All cuts were compen- 
sated for in the acceptance calculations. The effective masses for A hyperons 
were calculated by kinematic hypothesis only, by assigning a proton mass to 
the positive track and a pion mass to the negative track. 

Figure 4 shows the result of the effective mass calculation of the A ^ p + 7r~ 
hypothesis for the final selected data sample from the Au target. The dashed 
histogram represents the A effective mass reconstruction from the Monte-Carlo 
simulation of the experiment. The normalization per one central event is done 
for both data and Monte-Carlo. The small discrepancy in the tails is due to 
uncorrected distortions of tracks by the Au beam. The contribution of this 
discrepancy is well within the systematic errors. Different values in the peak 
of the distribution are due to different shapes of the momentum distribution 
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in the data and in the modeL Decay vertices with efTective masses in the 
range of 1.104-1.128 GeV/c^ were selected as A hyperons. The tails of the 
effective mass distributions were used for background subtractions. With the 
above mentioned cuts we obtained 3932 A hyperons above background for 
1.4 < y < 3.2 and rrit—mo < 0.7 GeV/(?. The A sample contains a contribution 
from higher mass states which decay to a A that passes the cuts. The main 
source of these is TP decay. This contribution is also contained in the model 
predictions. 




1.08 1.09 1.1 1.11 1.12 1.13 1.14 1.15 1.16 



M(7Tp) GeV/c' 

Fig. 4. Effective mass plot of the proton tt" hypothesis for 1.4 < y(A) < 3.2 and 
mt — mo < 0.7 GeV/c?. Data is shown by solid histogram, Monte-Carlo simulation, 
by dashed line. Both curves are normalized per one central event. 

The acceptance corrections of the distributions were calculated using a com- 
plete Monte-Carlo simulation of the effects of the apparatus and cuts on the 
final data sample. Events were generated using a modified HIJET model [13]. 
The (i/, mt) shape of the A spectrum in this model agrees with the measure- 
ments of experiment E810 [11] and the strangeness yield is normalized to K'^ 
and K~ measurements in Au + Au interactions [16]. The GEANT3 program 
was used to track the generated tracks through a magnetic field and the sim- 
ulated signals in the TPC were written out. The generated hits included all 
the known effects of detector apertures, efficiencies and distortions. (However, 
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see section 2.2 below). The results of this simulation were then analyzed by 
the same program used to analyze the actual data, including the tracking. A 
total of 25,000 central interactions were generated and analyzed. The accep- 
tance corrections on the final data sample were done as a function of rapidity 
and transverse mass. The A acceptance was in the range of 2-3% including a 
correction for neutral decays. 



2.2 Correction of Track Reconstruction Efficiencies in E891 



In the original analysis described above [10] the A hyperon candidates with 
protons having an azimuthal angle within the range, —52° < (p < 149°, were 
removed, where the azimuthal angle, </?, is measured with respect to y axis and 
the positive x direction corresponds to (f = 90°. The necessity of this angular 
cut was due to the fact that we could not simulate the dynamic component of 
the E and Ex B distortions caused by the positive ion cloud. The positive ion 
cloud charge density evolves on the time scale comparable to the duration of a 
spill (300 ms versus 1.2 s). The trigger time is randomly distributed relative to 
the beginning of the spill and the number of beam particles per spill fluctuates. 
This could not be successfully simulated in the Monte-Carlo. The angular cuts 
applied in the original publication did not guarantee that all the proton tracks 
will completely miss the distortion region. For these protons the Monte-Carlo 
tracks may be reconstructed with higher efficiency than the data. Since there 
is a correlation between the momentum of the A and the track coordinates in 
the TPC, this efficiency miscalibration can lead to a change of the measured 
A spectrum shape, which is crucial for determining the flow parameters. Using 
more stringent cuts to eliminate these tracks leads to a data set too hmited 
to determine the flow parameters. 

We have therefore recalculated the efficiency in the region of the distortion 
using a correction based on data from E866 [14,15]. The current analysis was 
also extended to lower rapidity than in [10]. The invariant differential spectra 
for p and tt"^ were measured in a wide range of rapidities and transverse 
momenta by experiment E866 and are found to be in a fairly good agreement 
with the model predictions by "A Relativistic Cascade Model" (ARC) [12]. 
We use the ARC model output for the positive particle spectra in order to 
determine the positive track reconstruction efficiencies as a function of the x 
and y coordinates of the track intercept with the front plane of the TPC: 

ATtracks f ■ I jyj events 



where A^^lrl number of positive tracks in the real data, N^^c;^^ that 
of the Monte-Carlo with ARC used as input and the numbers {i,j) indicate 
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Fig. 5. Positive track reconstruction efficiency as a function of the x and y co- 
ordinates of the track intercept with the front plane of the TPC. Beam impact 
coordinates are a; = 3.1 cm, y = 0.0 cm. 

the position of the 2 cm x 2 cm square to be caUbrated. The square size of 
2 cm is sufficiently small to insure a continuous dependence of efficiencies 
as a function of y and m^. We perform the cut on the impact parameter, 
bimp < 3 fm on the ARC event sample, to achieve the best description of 
our centrality selection criterion and to match the measured interaction cross 
section. Figure 5 shows the results of the efficiency calculation according to 
(1). One can see significant deterioration of the positive track reconstruction 
efficiency in the neighborhood of the positive ion sheet. The region plotted 
was chosen to best display the effect in the three dimensional plot. The fact 
that the efficiency goes smoothly to 100% away from the positive ion sheet 
demonstrates that the ARC model successfully extrapolated the E866 data in 
y and mj. 

Prom kinematics it follows that, for relativistic A hyperons, the momentum of 
a proton is almost collinear to the momentum of A which produced it. Thus 
the incident angle of protons from A's at the front plane of the first TPC 
is very close to that of a positive particle which came from the target. To 
illustrate the considerable overlap of the region of reduced efficiency of the 
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Fig. 6. Positive track reconstruction efficiency as a function of the x and y coordi- 
nates of the track intercept with the front plane of the TPC. The numbers in each 
square refer to the number of protons from A decay passing trhough that area. The 
contours represent the constant levels of efficiency ranging from 1.0, (outermost), 
to 0.4, (innermost), in steps of 0.075. 

positive track reconstruction with the area corresponding to the protons from 
A decay, we show the A counts (background subtracted) as a function of x, y 
coordinates in the front plane of the TPC in Figure 6. A significant fraction 
of the total number of A hyperons comes from protons passing through the 
region of reduced efficiency as can be seen from this figure. In the region of 
{x > —4 cm) and (y > —6 cm) the track efficiency drops below 80%. The 
number of A's going through this region, which have passed all cuts including 
the azimuthal angle cut, is about 1376. 

The calculation shows that A's affected by the inefficiency arising from the 
ion cloud have large rapidities and/or low transverse momenta (see Figure 7). 
The effective values of efficiency corrections for rapidity bins, integrated over 
mt, (1.7 < y < 2.0), (2.0 < y < 2.3), (2.3 < y < 2.6) and (2.6 < y < 2.9) 
are 1.03, 1.10, 1.21 and 1.36 respectively. The systematic errors for the A 
spectra, defined as the difference between the efficiency corrected data and 
the data without efficiency correction [10], do not exceed the statistical ones. 



10 



Track Efficiency Correction, E891 
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Fig. 7. A detection efficiency correction as a function of transverse kinetic energy 
and rapidity. 

The double differential spectra and rapidity distribution corrected for the track 
reconstruction efficiency are shown in Figures 9 and 10. 



3 The Blast Wave Pcirameterization 



An interesting phenomenon observed in relativistic Au + Au collision by ex- 
periment E891 is low-rrit suppression in the double differential spectra of A 
hyperons around midrapidity, a deviation from a single exponential scaling 
valid for lighter nuclei [11] and for pp interactions. The analogous effect has 
been observed by experiment E866 for protons [14]. It was suggested in [14] 
that the effect is due to the strong transverse expansion. The inffuence of col- 
lective nuclear flow has been estabhshed at different accelerator facihties. For 
instance, for the explanation of the different inverse slopes in the momentum 
distributions for protons and pions at Bevalac in the central Ne + NaF col- 
lisions at the lab energy of 800 MeV per nucleon, Siemens and Rasmussen 
[17] suggested that after thermal equilibration the nuclear matter experiences 
spherically symmetric expansion. According to the Bjorken scenario of one 
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dimensional longitudinal expansion [18], the spectra of particles produced in 
ultra-relativistic heavy- ion collisions are boost invariant. Such a spectrum can 
be represented as superposition of the thermal sources distributed uniformly 
over a limited interval of the longitudinal boost angle, r] [19] 



where rapidities y and t] are measured in the center of mass frame. The lim- 
its of the boost invariance interval, [—rimax,Vmax], are confined between the 
projectile and the target rapidities. Assuming that the longitudinal and the 
transverse expansion can be decoupled, P. Braun-Munzinger et al [19] used 
expression (2) integrated over the transverse mass squared to evaluate the 
rapidity distributions of different hadrons in Si + Al interactions at AGS en- 
ergy. Fairly good description of the rapidity distributions was achieved with 
the temperature, T = 120 MeV, and the mean longitudinal expansion veloc- 
ity, < vi >= 0.52. However, no attempt to parameterize the double differential 
distributions in a wide range of rapidities was made in this article. For the 
description of invariant differential multiplicity at midrapidity the authors ap- 
plied expression (A8) from [20], which is intended for the description of the 
transverse mass distribution integrated over rapidity, not the invariant differ- 
ential distribution! The disadvantage of such parameterization is that sources 
at different rapidities actually interfere with each other. Unlike the SPS ener- 
gies, where a modest boost of an observer in the longitudinal direction does 
not change the velocities of the receding nucleon pancakes, the AGS energy 
is not sufficiently high for the boost invariant expansion to take place. In ad- 
dition, nonzero curvatures of the reduced invariant differential multiplicities 
of baryons, S/rrit ■ d^N/dp^, on a logarithmic scale, indicate the presence of a 
transverse radial expansion with the strength varying as a function of rapidity. 
The longitudinal and the transverse expansion may be decoupled accurately 
only in the limit of nonrelativistic transverse velocities. Therefore, we believe 
that the correct procedure of determining the flow velocity profile is a pa- 
rameterization of the double differential spectra of particles. In this procedure 
the rapidity distributions and integrated transverse mass distributions will be 
described automatically. 

Usually the transverse velocity is extracted by fitting the spectrum assuming a 
power law of the transverse expansion velocity Vt oc and a constant density, 
p(rj) (see, for example [19]). From the cascade simulations it is known that 
the density resembles a gaussian profile, p{rt) = poc"'"'''^'^ , which can lead 
to significant consequences to the parameters obtained from a fit [22]. For 

^ The longitudinal boost angle, y = 0.5 In [(1 — Vz)], is essentially the ra- 

pidity of the moving source. 




rrit cosh (y—T)) /T 
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the traditionally used uniform density profile, as well as for a gaussian density 
distribution, the integration over the transverse expansion velocity selects non- 
zero velocities due to the rising factor, 2TiVt. The probability for a particle to 
have a transverse collective velocity, v^, has a peak at a value defined by 
parameters a and a in case of a gaussian density. Our data on the double 
differential A multiplicities, however indicates that a stronger selection of the 
transverse expansion velocities is present in reality, which is consistent with 
the creation of a transverse blast wave for the central Au + Au interactions. 
The development of the blast wave due to the sudden creation of hot dense 
matter was originally advocated in [17] where the spherical blast wave scenario 
was used to account for different apparent temperatures for p and tt. While 
apparent temperatures, the asymptotic values of the inverse slopes of reduced 
invariant differential cross sections, are not strongly sensitive to the probability 
distribution for collective velocities, the low-mt slopes are. It was pointed 
out in [23] that the step-like freeze-out distance distribution, widely used for 
experimental parameterization of transverse expansion (see for example [19]) 
is quite different from the one observed in RQMD, where the freeze-out close 
to the collision axis is suppressed. Before experiencing the final interaction the 
particles are propagated to the periphery of the system, where they acquire 
the essential collective motion [23]. This statement, together with the power- 
like transverse velocity profile is consistent with a blast wave scenario. It was 
also claimed in [23] that it is impossible to describe a spectrum with a single 
temperature and transverse expansion velocity. We will demonstrate that the 
spectrum of A hyperons can be described in terms of an aspherical blast wave 
with a single temperature. 

The influence of the transverse blast wave on the low-m^ slopes of the reduced 
invariant differential cross sections can be understood by noting that the spec- 
trum in a given rapidity bin originates primarily from sources, macroscopic 
elements of nuclear matter, with a limited rapidity range. One can thus fol- 
low this influence by considering the transversely expanding ring at this given 
rapidity. Boosting the invariant Boltzmann distribution in the transverse di- 
rection if by velocity 

^Oiv^ = _ COS <^)e-^(^-f*--^)/^, (3) 



and integrating this expression over the azimuthal angle, (/?, one can derive 
the equation for a spectrum from the transversely expanding ring: 



= 27r7{^:/o(0 -P*^^i(0}e- - , (4) 
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where 



C = IPtv/T, 7 = l/Vl - v^, (5) 



2-ir 27r 

and functions Io{x) = (27r)-i / e'=''°^'^dip and 7i(x) = (27r)-^ / cos(^e^'=°^'^d(/? 



are the modified Bessel functions of the zeroth and the first order respectively. 
The low-mt slopes for the of reduced invariant differential cross section can be 
found by resolving the logarithm of equation (4) divided by rrit into a Taylor 
series in the vicinity of — >• mo, assuming S = rrit, and saving only the linear 
term: 

, f S d^N\ -fi'ymov^ - 2T(1 + v^)} , , 
'"Urfi? 2T^ ^K-m.). (6) 
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Fig. 8. The low-mj slopes for various hadrons as a function of the transverse expan- 
sion velocity for a fixed temperature of T = 100 MeV. 

The factor multiplying the transverse kinetic energy, rrit — ''^Oi represents the 
low-mt slope of the reduced invariant differential cross section on a logarith- 
mic scale. On Figure 8 the low-m^ slopes are shown for various hadrons as 
a function of the transverse expansion velocity for a fixed temperature of 
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T = 100 MeV. One can see that for sufficiently high transverse expansion 
velocities the slopes become positive, yielding the low-m^ suppression. It is 
worth noting that for tt mesons the low-mj slope shows the opposite behav- 
ior, though of smaller magnitude. Generally, the condition mo > 2T must be 
satisfied for low-m^ suppression to take place. The higher the temperature of 
the hadron gas the less affected is the low-mt slope of the spectrum. 

We represent the inclusive spectrum as a continuous superposition of flowing 
Boltzmann sources with an effective temperature, T. The effective flow is 
described by an azimuthally symmetric surface in the space of the longitudinal 
boost angle, 77, and the transverse two-dimensional velocity, the Blast Wave. 
The blast wave in our approach belongs to the interval, [—f]max,Vmax], where 
both limits are measured in Au + Au center of mass system. An expression 
for the invariant differential multiplicity reads as follows: 



d^N r 

^-^= J dr]Sf{y,r],mt), (7) 

Vmax 

SfiV: V: m) = -A{i^){xm) - (8) 

TT 

^ 7^(77) mt cosh {y - ?]) ^ 7t(^) Pt Vtjv) 



Vmax 



Choosing the flow velocity proflle and the rapidity distribution of the emission 
power of the blast wave, we employ the following simple functions: 

vtiv) = vo{i - vVvL.r, Mv) = Ao{i - vVvL.r, (11) 



where a, u > 0. Rapidities, y and rj, arc measured in the center of mass frame 
of the collision system and the transverse velocity, Vt, is measured in a system 
co-moving with 77. The functions Vt{r)) and A{'r]) are chosen to allow one to 
tune the shape of the blast wave to the data. One can easily see that with 
z/ = 1/2 and Vq = r]max the function Vt{ri) describes a spherical blast wave, 
where for a = the emission power is uniformly distributed along the surface 
of the blast wave. With rjmax ~^ and Vq > 0, one encounters purely transverse 
expansion, as in the Landau scenario [21]. It can be also seen that the boost 
invariant Bjorken cylinder can be described with 1/ — a — 0. With positive 
values of parameters u and a, the functions above reflect the experimental fact 
of increase of the transverse expansion velocity towards midrapidity, which is 
seen from the increase of the low-m^ slope of the reduced invariant differential 
cross section on a logarithmic plot (see Figure 9). 
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A Spectra Au+Au, E891 




Fig. 9. Double differential multiplicities for A production in central Au + Au inter- 
actions measured by experiment E891. The models shown are a) RQMD cascade, 
by dashed lines, b) RQMD with mean fields, by dotted lines c) Blast Wave fit, by 
solid lines. 



Mean rapidity and the transverse velocity in the blast wave are given by the 
following expressions: 



1 

Jx{l-x^)-dx T{a + 3/2) 





1 



^ r(a + z/ + l) r(a + 3/2) 

kl-x^)-dx r(a + ^ + 3/2) r(a+l) 



where averaging is performed over the forward hemisphere in the center of mass 
system. The rapidity distribution and the total number of particles emitted 
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Au+Au - A Rapidity Distribution, E891 
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Fig. 10. Rapidity distribution for the A hyperon production in central Au + Au in- 
teractions measured by experiment E891. The models shown are a) RQMD cascade, 
by dashed lines, b) RQMD with mean fields, by dotted lines c) Blast Wave fit, by 
solid lines. 

by this source are 



dN{y) mle t 
dy 



rriQ Tjmax 



TT J 



dNboitzjy - rj) 
dy 



(14) 



Aom^rjmaxe— T{a + 1) 
^ T{a + 3/2) 



holtz -I 



(15) 



respectively, where V{x) — \ *dt, (x > 1), is the gamma function, and 





2 • T 



dy 



mo cosh (y) ml cosh (y) 



ITIQ cosh(i;) 

e T (^16) 



is the Boltzmann rapidity distribution. 
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RQMD mean fields 


{(jev /c ) 
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4.0D It U.OZ 




U.OU It U.Uo 


U.o4 It U.UD 


u.4y It U.Uo 


rv 

(_x 


n c:r)_|_1.17 
U.OUitQ 5Q 


n qo_|_1.4 
u.yoitQ gg 


n 99_|_0.31 
U.ZZitQ 22 


V 


0.75ibHI 


0.98ib3,-^§ 


0.57 lb 0.22 


T (MeV/c^) 


96 lb 37 


136 lb 33 


95 lb 13 


^TTLdX 


1.36 lb 0.38 


1.1 ±0.30 


1.18 lb 0.10 


y'^/NDF 


32/43 


58/58 


60/47 


<Vt> 


0.40 lb 0.12 


0.27 lb 0.08 


0.40 lb 0.04 


<rj> 


0.58 lb 0.07 


0.42 lb 0.02 


0.54 lb 0.02 


Na 


20 lb 3.33 


18 lb 1.52 


19 lb 1.12 



Table 1 

Results of the Blast Wave parameterization of A double differential multiplicities 
for E891 data and RQMD predictions. Six main parameters are listed first. The 
mean transverse expansion velocity, < vt >■, longitudinal boost angle, < r] >, and 
mean multiplicity, A^a, are obtained using formulas (12), (13) and (15). 

We have applied the formulas listed in this paragraph to obtain the expansion 
parameters as well as the mean temperature of A hyperons after the freeze-out. 
The obtained parameters Aq, vq, a, z/, T, rjmax are listed in Table 1 together 
with the mean source rapidity, < rj >, the mean transverse expansion velocity, 
< Vt >, and the integrated yield, A^a derived from the above six parameters. 



4 Results and Discussions 

In Table 1 the results of the parameterizations of the double differential mul- 
tiplicities for the A production measured by E891 are listed together with the 
results obtained from the cascade version of RQMD and the version of RQMD 
with mean field potentials turned on. For comparison with RQMD model pre- 
dictions we have selected the sample of central events with the cut on the 
impact parameter, 6 < 3 fm. This value of the cut was chosen, since the re- 
sulting Au + Au geometric cross section, vr6^^^, matches the cross section we 
obtain and the resulting vr" differential multiplicities agree well with our mea- 
surements 0. One can see that the cascade version of RQMD underpredicts 
the longitudinal and the transverse flow, which can be verified by comparing 

^ In Reference [10] the cut of 6 < 4 fm was chosen in models for comparison with 
our data. 
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the mean transverse expansion velocity, < Vt >, and the mean longitudinal 
boost angle, < 7] >. This fact can also be seen on the plots of double dif- 
ferential multiplicity for A production (Figure 9) where the cascade version 
of RQMD fails to describe the low-mj suppression effect in the midrapidity 
region. The rapidity distribution of the A hyperon measured by E891 (Figure 
10) has a larger width than predicted by RQMD cascade, which corresponds 
to the smaller strength of longitudinal flow in the cascade. One can also see 
in Table 1 that all the parameters of the Blast Wave obtained by fitting the 
results of RQMD with mean fields turned on agree with the ones obtained 
by fitting the data within experimental uncertainties. The low-m^ suppression 
effect is reproduced by this version of RQMD and the width of the rapidity 
distribution agrees well with the experimentally measured width. 



5 Conclusions 

Based on the agreement of the predictions of RQMD with mean fields with 
our measurements for the double differential and rapidity distributions for A 
production in the central Au + Au interactions and underestimation of the 
fiow strengths by RQMD cascade we can confirm that the mean field effects 
play important role in the dynamics of the Au + Au interactions at the AGS 
energy. 
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